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Preface 


This  paper  applies  the  simple  Taylor  series  ap¬ 
proximation  to  the  intricate  problem  of  correcting 
inaccurate  probe-positioning  in  planar,  near-field  scanning. 

This  is  not  to  imply  that  the  implementation  is  as 
simple  as  the  concept.  Throughout  the  past  eight  months  1 
have  learned  that  simple  concepts  can  be  quite  difficult  to 
implement . 

There  are  many  subtleties  that  the  student  new  to  this 
area  encounters.  Dealing  with  the  minor  differences  between 
todays  modern  main-frame  computers  alone  can  delay 
indefinitely  the  more  improtant  development  and 
interpretation  of  the  computer  analysis. 

There  are  several  individuals  who  helped  me  overcome 
these  differences.  Particular  thanks  goes  to  my  sponsor, 

Dr.  Arthur  D.  Yaghjian,  RADC/EEC,  Hanscom  AFB.  Without 
his  help  this  thesis  would  not  have  been  possible.  Thanks 
also  goes  to  my  thesis  advisor,  ILt  Randy  Jost,  for  his 
guidance  and  patience  to  the  point  of  correcting  my 
spelling  mistakes. 

Finally,  I  would  like  to  acknowledge  Allen  C.  Newell, 
of  the  National  Bureau  of  Standards,  as  the  one  who  first 
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proposed  this  method  of  probe  position  compensation. 
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Abstract 

In  recent  years,  near-field  antenna  measurement 
techniques  have  gained  a  good  deal  of  acceptance.  There 
are  several  errors  in  the  computed  far-field  patterns  caused 
by  measurement  inaccuracies  in  the  near-field  data.  This 
paper  deals  with  deterministic  errors  introduced  by  probe¬ 
positioning  errors  in  planar,  near-field  scanning. 

By  utilizing  basic  near-field  theory,  as  well  as  a 
knowledge  of  the  positioning  errors,  it  is  possible  to 
estimate  the  fields  at  the  correct  probe  position.  A 
computer  program  which  lessens  the  effect  of  probe 
positioning  errors  by  means  of  a  truncated  Taylor  series 
expansion  is  used  to  demonstrate  this  improvement. 

To  simulate  inaccurate  probe  positioning,  a  database  of 
calculated  near-field  values  for  a  linearly  polarized, 
uniformly  illuminated,  circular  aperture  was  used.  Next, 
the  position  correction  program  was  applied  to  computer 
generated  inaccurate  data.  Finally,  far-field  patterns 
were  calculated,  and  compared  using  both  the  corrected  and 


uncorrected  data.  Results  of  the  comparisons  are  presented. 
Limitations,  and  areas  of  application  of  this  routine  are 
discussed. 
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Near-field  scanning  techniques  are  quickly  becoming 
accepted  as  an  efficient,  and  accurate  method  for  the 
determination  of  antenna  patterns  [9:492].  Near-field 
techniques  offer  many  advantages  over  conventional  antenna 
measurement  facilities  including:  all  weather  operation, 
reduced  ambient  interference,  and  security  for  delicate,  or 
classified  apparatus  [11:101].  To  further  enhance  the 
reliability  of  antenna  parameters  predicted  by  this 
procedure,  any  source  of  error  must  be  acknowledged  and 
limited  where  possible. 

Near-field  scanning  requires  the  measurement  of  probe 
output,  in  both  phase  and  magnitude,  over  a  predetermined 
scanning  surface  such  as  a  sphere,  circular  cylinder,  or 
plane  [7].  This  is  accomplished  by  positioning  a  probe  at  a 
point  on  the  scanning  surface  through  the  use  of  mechanical 
devices.  Since  mechanical  devices,  such  as  positioners, 
are  not  infinitely  accurate,  some  error  in  the  probe 
positioning  exists.  Inaccurate  probe  position  specification 
leads  to  erroneous  far-field  patterns,  calculated  from  the 
incorrect  data  gathered  by  the  probe. 

Problem  Statement 

This  thesis  will  examine  sensitivity  to  inaccurate 
probe  positioning,  as  well  as  a  possible  method  for  reducing 


that  sensitivity  by  means  of  a  truncated  Taylor  series 
approximation  of  the  near-field  data. 

Scope 

A  computer  program*  created  to  generate  the  near-field, 
transverse  electromagnetic  fields  on  a  plane  in  front  of  the 
antenna  is  used  as  a  data  base  for  the  computation  of  far- 
field  antenna  patterns,  inaccuracies  in  probe  positioning 
are  simulated  by  the  manipulation  of  this  program.  This 
data  set  is  used  to  calculate  the  far-field  pattern  of  the 
test  antenna  using  planar  scanning  techniques  explained  in 
Chapter  11.  The  resulting  pattern  is  then  compared  to  the 
pattern  from  a  perfect  scan  for  various  degrees  of 
inaccurate  positioning. 

Also,  the  position  correction  routine  is  used  to 
compensate  for  erroneous  probe  positioning  and  again  the 
far-field  pattern  is  calculated.  This  "corrected”  pattern 
is  compared  to  the  uncorrected  pattern  as  well  as  the 
pattern  from  the  perfect  scan,  to  determine  whether  or  not 
any  improvement  has  occured.  This  procedure  is  also 
repeated  for  various  degrees  of  positioning  error. 
Assumptions 

To  lessen  the  number  of  the  computations,  and  to 
conserve  computer  time,  two  assumptions  are  made: 

1)  The  simulated  positioning  errors  are  discrete  in 
nature,  and  consist  of  a  single  displacement  at  the  center 
line  of  the  scanning  surface. 


2)  The  hypothetical  probe  used  is  an  ideal  electric 
dipole  sampling  the  electric  fields. 


This  research  effort  began  with  a  review  of  the 
development  and  theoretical  formulation  of  planar,  near¬ 
field  scanning  techniques.  The  results  of  this  review  that 
are  pertinent  to  this  presentation  are  given  in  Chapter  II. 

After  an  understanding  of  the  basic  theory  of  planar, 
near-field  scanning  was  obtained,  a  data  base  of  near-field 
values  was  needed  to  investigate  the  effect  of  probe 
positioning  errors.  A  computer  program  that  calculates  the 
fields  of  a  circular  aperture,  linearly  polarized  antenna  as 
a  function  of  frequency,  aperture  radius,  and  distance  from 
the  antenna  was  used  to  generate  this  data.  A  description 
of  that  program  is  presented  in  Chapter  HI. 

Next,  a  program  which  calculates  the  far-field  pattern 
of  an  antenna  from  its  near-field  values  was  created  to 
examine  the  effect  of  inaccuarate  probe  positioning.  The 
position  correction  routine  is  an  integral  portion  of  this 
program.  Both  procedures  are  explained  in  detail  in 
Chapter  IV. 

Finally,  the  effect  of  positioning  errors,  as  well  as 
the  effectiveness  of  the  position  eorrection  routine  in 
limiting  these  errors  is  examined  through  a  computer 
analysis  of  the  far-field  patterns  produced.  This  analysis, 
and  the  results  it  provided  are  presented  in  Chapter  V. 


In  this  chapter,  the  basic  theory  behind  planar,  near¬ 
field  measurements  is  developed.  The  definitive  analysis  in 
this  area  was  performed  by  O.  M.  Kerns  [12].  The  approach 
used  involves  the  expansion  of  a  set  of  measured  field 
values  into  a  summation  of  its  elementary,  planar  wave 
f unct ions . 

The  following  is  a  list  of  the  variables  and  the 
notations  that  will  be  used  throughout  the  theoretical 
development  of  this  report. 

£(r)  =  The  electric  field,  magnitude  and  phase  of  each 

component,  as  a  function  of  position. 

A(k)  =  The  vector  amplitude,  magnitude  and  phase  of  each 
component,  as  a  function  of  propagation  direction. 

7"  =  The  position  vector  x 'a'jj  +  y  'a'y  +  z 

k  =  The  propagation  vector  kjj  ^  +  ky  'a'y  +  k^  'a'^* 

=  The  frequency  in  rad/sec. 

c  =  The  permittivity  of  the  propagation  medium. 

=  The  permeability  of  the  propagation  medium. 

R  =  The  tangential  position  vector  x  ^  y 

"k  =  The  tangential  direction  vector  kj^ 'a'jj  +  ky  ^a^. 

dK  =  The  two  dimensional  elemental  patch  in 
k-space,dkjjdky. 

dR  =  The  two  dimensional  elemental  patch  in 
real  space,  dxdy. 

k  =  The  phase  propagation  constant  In/K  . 


The  propagation  of  electromagnetic  waves  within  free- 
space  is  described  by  Maxwell's  equations.  These  equations 


lead  directly  to  the  vector  Helmholtz  equation  [13:124]. 

v2e  (F)  +  k2  e  (F)  *  0  (1) 

The  simplest,  non-trivial  solution  to  equation  (1)  for 
a  homogeneous,  isotropic,  source  free  region  (z  >  0)  is 
given  by  [1:145]: 

E  (r)  =  A  (F)  e~^  (2) 


where:  T*  is  the  position  vector 

IT  is  the  direction  of  propagation 

A(k)  is  the  vector  amplitude 

of  the  wave  as  a  function  of  the 
direction  of  propagation. 

By  substituting  equation  (2)  into  equation  (1),  we  obtain 
[1:144] : 


+  kj  +  k 


.2  =  k2 


(3) 


An  interesting  property  of  equation  (3)  is  that  given  two 
components  of  the  propagation  vector,  say  k,^  and  ky,  the 
third  component,  k2i  is  specified. 

+  ky^  +  kz^  *  (4) 

All  values  of  k^,  and  ky  which  yield  a  real  value  for  k2 
define  the  propagating  modes.  However,  modes  for  which 
equation  (4)  produces  a  complex  result  are  evanescent,  or 


non-propagating.  Still  these  evanescent  modes  must  satisfy 
the  convergence  condition»  namely  [5:374]: 


lim  A  (k)  e-  i  a  0  (1 

f ‘►OP 

It  is  for  this  reason  that  we  choose  the  following  [1:146; 


10:499] : 


'J  (k^  -  K^)  k2  ‘i  k2 

-i  y  (k^  -  k2)  k2  >  k^ 


Note  that  k^;  is  either  positive  real,  or  negative  imaginary. 

Since  we  have  assumed  a  source  free  region,  Maxwell's 
equations  require  that  [13:112]: 


V  •  E  (F)  =  0 


substituting  equation  (2)  into  equation  (7)  yields: 

]r  •  A(k)  =  0 


Equation  (8)  implies  that  given  the  direction  of 
propagation,  and  two  components  of  the  vector  amplitude,  the 
third  component  is  uniquely  defined  as: 


*  ^y'^y^ 

For  any  complex  system,  the  complete  electric  field,  at 
any  point  in  the  source  free  region,  can  be  expressed  as  a 
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summation  of  the  elementary  wave  functions  over  all  possible 
propagation  directions  [5:374;  15]. 

E(r)  ■  2ir  A(]c)e“^  3k"  (10) 

If  A(k)  in  equation  (10)  is  expressed  in  rectangular 
components,  the  integral  can  be  reduced  to  three  scalar 
expresions,  corresponding  to  each  of  the  unit  vector 
components. 

(11) 

A  closer  examination  of  the  integral  expressions  given 
in  equation  (11)  reveals  that  they  are  the  two  dimensional 
Fourier  transforms  of  the  components  of  X(k)  [4:381].  The 
inverse  relation  is  given  by  [5:375;  15]: 

m  m 

A(lt)  -  -Ly"  f  ^<ir  (12) 

2w 

By  specifying  the  xand  y  components  of  the  electric 
field  over  a  plane  defined  by  z=Zq  we  can  determine  from 
equation  12  the  x  and  y  components  of  A(k),  and  by  using 
equation  (9)  we  can  uniquely  specify  the  remaining 
component  A-* 


(13) 


Given  "ACk),  the  electric  field  is  known  for  all  space 
by  equation  (10).  Generally,  this  equation  must  be 
performed  on  a  computer.  However,  for  the  case  of  the  far- 
field  pattern,  equation  (10)  can  be  evaluated  analytically 
using  the  method  of  steepest  descent  [5:375;  3:1675;  6:253]. 


Epp( r ) 


1  2ii  kio  A(ko)  - 


(14) 


where:  l^o  ~  ^  ^r 

At  =  sin(^  )cos(0)x‘  +  sin(d  )sin(0)^  +  cos(0  Yz 
Finally,  note  that  since  kjjo^  +  =  ic2  sin2(e  )  <  k^ 

equation  (14)  involves  only  the  spectrum  of  the  propagating 
modes . 

CLabfiJBa&iiiaaiaK_££.Laj:& 

The  foregoing  development  assumed  that  the  tangential 
field  was  known  over  a  plane  defined  by  z=Zq.  It  is 
impossible  to  create  a  perfect  measuring  apparatus, 
therefore  some  error  is  always  present  in  the  form  of 
inaccurate  probe  positioning.  This  error  can  occur  in  two 
distinct  forms: 


1.  Random  error  introduced  by  loose  machinery, 
temperature  variations,  or  ground  vibrations. 

2.  Deterministic  positioning  errors  due  to  limited 
machine  accuracy. 

The  first  is  time  variant,  and  can  be  eliminated  by  making  a 
statistical  average  of  several  scans.  For  a  given  system, 
the  deterministic  errors  can  be  measured  with  greater 
accuracy  than  the  positioning  machinery  is  capable  of 
achieving,  and  it  is  this  type  of  positioning  error  that  is 
the  topic  of  this  investigation. 

The  deviation  of  the  probe  from  the  theoretical  plane 
of  measurement  will  introduce  errors  in  the  far-field 
patterns.  To  counter  this  we  will  use  the  spatial  deviation 
measured  and  a  knowledge  of  the  dependence  of  the  tangential 
fields  on  variations  in  the  z  direction  to  estimate  the 
value  of  the  fields  at  the  theoretical  plane  of  measurement. 
We  will  consider  errors  in  the  z  direction  only. 

To  estimate  the  value  of  the  tangential  electric 
fields,  we  will  make  use  of  a  truncated  Taylor  series 
expansion.  The  general  form  of  the  Taylor  series  is  given 
by  [18,  305]: 


f(z) 


y'  ^  d  ^ ^ 

nao  dz^"J 


f(.)| 

z«zo 


(15) 


If  we  assume  that  the  higher  orders  of  the  series  are 
dominated  by  the  first  and  second  terms  we  can  say: 


E(z)  •  ^(zo)  4-  I  E(2)|az 

2»Z0 

By  examining  equation  (10)  we  can  see  that: 


(16) 


^  an  A  He)  e  (17) 

To  =  xajj  +  yay  +  Zoa^ 

So,  if  we  knew  the  vector  amplitude  function  of  a  perfect 
scan,  A(k),  we  could  determine  the  first  partial  of  the 
electric  field  with  respect  to  z  and  approximate  the  value 
of  the  fields  on  the  plane  z=Zq.  A(k)  is  not  yet  known,  but 
we  can  approximate  A(k)  by  using  equation  (13)  to  calculate 
a  ^(k)  which  is  the  vector  amplitude  function  calculated 
using  the  measured  data,  and  neglecting  the  positioning 
errors  . 


Ao(ir) 


2. 

2r 


(Zl )® 


+i 


(18) 


Therefore,  our  equation  for  the  tangential  fields  on  the 
theoretical  plane  of  measurement  becomes  [17:24]: 


(zo) 


(zl ) 


■  T  /  J^z)  Ao 


2w  -• 


(k)e  “  (19) 


In  general,  4z  =  (zj  -  Zq)  is  dependent  on  the  x,  y 
coordinates  of  the  point  in  question. 


Before  a  study  into  probe  positioning  errors  can  begin, 
a'database  of  near-field  electric  field  values  along  with 
their  accompanying  far-field  patterns  must  be  obtained.  For 
this  investigation,  the  database  must  include  flexibility  of 
the  positioning  discrepancies  as  well  as  a  reliable  degree 
of  accuracy  so  as  not  to  mask  the  outcome  of  any  correcting 
routines  used.  Finally,  since  the  database  will  be  accessed 
repeatedly,  it  should  be  as  easy  to  use  as  possible. 

For  these  reasons,  a  computer  program  that  calculates 
the  tangential  electric  fields  in  the  near-field  of  a 
linearly  polarized,  uniformly  illuminated,  circular  aperture 
antenna  was  chosen  as  a  source  of  raw  data  over  actual 
measurments.  The  original  analysis  and  program  was  written 
by  Dr.  Arthur  Yaghjian  RADC/EEC,  Hanscom  AFB,  and 
provides  data  based  on  parameters  supplied  by  the  user. 

The  main  parameters  are,  ka  (i.e.  27ra/A.)  and  z/a  where  z 
is  the  distance  to  the  plane  of  computation. 


General  Outline 

To  create  the  database  program,  the  fields  of  the  test 
antenna  were  expanded  into  elementary  spherical  harmonics 
[8].  This  was  accomplished  by  first  examining  the 
corresponding  acoustical  problem  of  the  two-sided  piston 
radiator.  Since  the  differential  equations  of  acoustics  are 
only  slightly  different  from  those  governing 
electromagnetics,  similar  radiating  systems  yield  similar 
results.  By  expanding  the  acoustical  pressure  pattern  of 
the  two-sided  piston  radiator  into  spherical  harmonics,  an 
expression  is  obtained  that  can  be  used  as  a  starting  point 
for  the  more  difficult  electromagnetics  problem.  This 
process  leads  to  an  equation  that  is  an  infinite  sum  of 
individual  spherical  modes,  an  equation  that  lends  itself  to 
computation  by  computer. 

The  geometry  that  the  database  program  is  based  on  is 
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The  equation  which  yields  the  electric  field  values  for 
all  points  outside  the  source  sphere  (i.e.  |T|  >  a)  is  given 


by  [8:13]: 


Ex(r) 


-(ka)2 

ITTcr 


•.io 

even 


t  (i+ixcose  K; 

+  Cose  (till  -  hji+i)PA  1 


)  YY  ('^"1-2,8)  1  I  (t-H-2a)  i  1  (-1  ) 

2**’  q-l\  -q)  t  +q)  i  J 

X  (jq  (l«j  j-q  (j^)J 


(A-I)ll 

(A)  1 
”2 


(-1) 


£ 

2 


(20) 


Note  that  ail  of  the  summation  constants,  ,  are  real. 
Therefore  the  phase  information  is  given  by  expressing  the 
Hankel  functions  as: 

^n^=  in  ~  *  "n  <21) 

At  first  glance  it  would  appear  that  this  equation 

would  supply  both  the  near>field  values,  and  the  far-field 
values.  Practically,  this  is  not  the  case  due  to  the  fact 
that  as  the  computation  surface  is  moved  away  from  the  test 
antenna,  the  number  of  terms  needed  to  insure  accuracy 
becomes  prohibitivly  large.  Fortunately,  this  is  not  a 


problem  since  the  far-field  pattern  of  this  system  can  be 
expressed  in  closed  form.  The  far-field  pattern  is  given  by 


[2:417]: 
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Near-Field  Computation 

The  near-field  database  program  addresses  each 
component  of  equation  (20)  separately.  There  are,  however, 
four  special  functions  that  are  common  through  out  the 
computational  procedure  of  a  single  data  set.  These  special 
functions  are  calculated  in  subroutines  called  by  the  main 
program  and  they  include: 


1)  The  Factorial  Function 
Calculated  in  the  subroutine  Fact 

2)  The  Fac-Factorial  Function 
Calculated  in  the  subroutine  Facfact 

3)  The  Spherical  Bessel  Functions 
Calculated  in  the  subroutine  Sphebe 

4)  The  Legendre  Polynomial  Functions 
Calculated  in  the  subroutine  Legend. 


Each  subroutine  is 
that  define  the  special 


governed  by  recurrence  relations 
function  that  it  calculates. 


The  Subroutine  Paet 


This  subroutine  is  a 


simple  algorithm  that  calculates 


the  factorials  of  the  first  N  integers  where  N  is  a  limiting 
value  supplied  by  the  main  program.  The  equation  used  to 


calculate  these  values  is  given  by: 


P 

P!  -  77  n  (0!=  1)  (23) 

n«l 

The  Subroutine  Faefaet 

This  routine  calculated  the  double  factorials  of  the 
first  2N-I-1  integers  where  N  is  a  maximum  value  supplied  by 
the  main  program,  and  is  based  on  the  parameters  supplied  by 
the  user.  The  expression  used  to  calculate  these  values  is 
given  by: 

P 

(2P+1)1!  *  77  (2n+1)  . 

n-l  <24) 


The  Subroutine  Sphebe 


This  routine  is  slightly  more  complex  than  the 
preceding  two,  since  it  must  calculate  the  spherical  Bessel 
functions  of  an  arbitrary  order  and  arbitrary  arguments. 
That  is,  it  must  calculate  the  ordinary  spherical  Bessel 
functions  as  well  as  the  spherical  Neumann  functions.  The 
equat ion  used  is : 


(x)  +  fn+i  (*)  - 


(25) 


All  cases  use  known  forms  of  jn(x)  and  n^(x)  or  the 
above  recurrence  relation. 

For  j„(x)  we  use  a  descending  computation  starting  at 
N>n  and  set  JnC*)  *  jN-l<*^  -  ldl  <<  1 

For  n|,(x)  we  use  an  ascending  computation  using  the 
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known  forms  of  the  first  two  Neuman  functions. 
The  Subroutine  Leyend 


This  routine  calculates  the  Legendre  Polynomials  up  to 
order  N,  and  for  argument  cos(^)t  where  N,  and  cos(6)  are 
provided  by  the  main  program.  The  algorithm  is  based  on  the 
equat i ons : 

(2n+i)  X  P(n)(*)  -  (n+1)  P(n+1)(x)  +  n  P(n-1)<*) 


Pq  (Cos0)  =  1  ^26) 

P]  (CosB)  *  Cose 
The  Main  Program 

The  main  program  begins  calculation  of  field  values  by 
first  calculating  the  summation  constants  given  in 
equation  (20).  Next,  the  program  calculates  the  phase  and 
magnitude  of  the  fields.  This  is  done  by  taking  full 
advantage  of  the  radial  symmetry  of  the  geometry.  Examining 
equation  (22),  it  is  noted  that  the  far-field  pattern  is 
symmetric  about  the  z-axis.  Therefore,  the  program  only 
needs  to  calculate  the  neer>field  values  on  any  ray  normal 
to  the  z-axis,  and  lying  in  the  plane  of  computation.  This 
reduces  the  number  of  values  to  be  calculated  and  the  time 
of  computation. 

The  main  program  introduces  a  variable  R,  which  is  the 
distance  off  of  the  bore  sight  of  the  antenna  on  the  plane 
of  computation.  It  then  calculates  the  phase  and  magnitude 
of  the  electric  field  at  that  point.  The  equation  for 


these  values  is  given  by  [8:15]: 
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Finally,  the  program  writes  these  values  to  a  data  file 
which  is  used  as  the  input  data  to  the  Far-Field 
Transformation  Program. 


IV. 


Once  a  database  of  unprocessed  values  has  been  obtained 
for  the  test  antenna,  the  next  step  is  to  create  programs 
which  reduce  this  data  into  a  form  suitable  for  comparison. 
This  is  accomplished  through  Transformation  Routines. 

These  programs  use  the  near-field  values  calculated  in 
the  database  program  to  compute  the  vector  amplitude 
function,  A(k)  through  a  numerical  evaulation  of  equation 
(13).  This  data  is  then  used  to  evaluate  equation  (17) 
which  yields  the  first  partial  derivative  of  the  electric 
field  with  respect  to  tne  variable  z,  at  the  plane  of 
measurement.  Next,  the  original  near-field  data  is  combined 
with  its  first  partial  by  implementing  equation  (19)  and  the 
vector  amplitude  function  is  recalculated  using  this 
corrected  data  set.  Finally,  the  far-field  pattern  is 
plotted  as  a  function  of  $f  and  compared  to  the  pattern 
calculated  from  a  perfect  scan.  Figure  2  is  a  block  diagram 
of  this  procedure. 


tr 


wf 
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The  process  depicted  in  figure  2  is  performed  in  four 
steps.  This  is  necessary  due  to  hardware  limitations  on 
the  maximum  job  size  presented  by  the  main  frame  computer 


used.  Each  step  represents  a  separate  computer  program,  and 
each  program  evaluates  one  of  the  equations  outlined  in 
Figure  2.  The  programs,  in  order  of  operation,  are: 

1)  The  Transformation  Program,  TRPROG 

2)  The  Derivative  Program,  DRPROG 

3)  The  Combination  Program,  OOMPROG 

4)  The  Output  Program,  OUTPROG 

The  following  is  a  detailed  description  of  each  program. 
TRPROG 

This  program  evaluates  equation  (13)  by  first  reading 
in  the  near-field  data,  and  then  performing  the  double 
integral  for  each  point  in  the  transform  space  (K-space) 
where  a  value  for  the  vector  amplitude  function  is  required. 

There  is  a  complication  in  reading  the  near-field 
values  supplied  by  the  database.  As  described  in  Chapter 
II,  the  database  program  calculates  the  near-field  values 
for  a  linearly  polarized,  uniformally  illuminated,  circular 
aperture  antenna.  It  accomplishes  this  taking  full 
advantage  of  the  cylindrical  symmetry  of  the  problem. 
Consequently,  the  final  output  of  the  program  is  in 
cylindrical  coordinates.  Before  this  data  can  be  used  by 
the  TRPROG,  it  must  be  expressed  in  cartesian  coordinates. 
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To  acheive  this,  the  program  first  creates  a  two  dimensional 
array  into  which  the  data  is  to  be  read.  This  array 
represents  a  scanning  plane,  and  each  storage  location  has 
associated  with  it,  a  pair  of  cartesian  coordinates 
(Xi,  y|).  Next,  TRPROG  reads  the  cylindrical  data  from  the 
database  and  stores  it  in  two  one-dimensional  arrays,  one 
for  the  phase,  and  one  for  the  magnitude.  Each  storage 
location  of  the  one  dimensional  arrays  has  a  radial 
coordinate  associated  with  it.  Since  the  electric  field  of 
the  test  antenna  is  circularly  symmetric  about  the  z-axis, 
only  the  radial  component  is  required  to  specify  the  field 
values  for  a  given  location.  Finally,  using  the  cartesian 
coordinates  associated  with  each  storage  location,  the 
radial  distance  is  calculated.  The  phase  and  magnitude 
values  from  the  database  whose  associated  radial  dimension 
most  closely  matches  this  calculated  radial  dimension  is 
read  into  the  storage  location  as  a  complex  value,  see 
figure  3. 


Figure  3.  The  Coordinate  Transformation  Process. 


After  the  near-field  data  has  been  read,  TRPROG 
perforins  the  double  integration  of  equation  (13).  This 
entails  a  double  integral  for  each  value  of  the  vector 
amplitude  function  to  be  calculated.  If  one  examines 
equation  (17),  the  next  equation  to  be  implemented,  it  is 
evident  that  a  plane  of  vector  amplitude  values  will  be 
required  for  the  next  computation.  This  determines  which 
vector  amplitude  values  must  be  calculated  in  TRPROG.  The 
simplest,  and  most  direct  method  of  computation  is  to  create 
a  two  dimensional  array  similar  to  that  containing  the  near¬ 
field  data  and  associate  each  storage  location  with  a  pair 
of  coordinates  from  K-space. 

Computation  time  is  a  major  concern  in  choosing  a 
method  for  calculating  the  double  integrals  of  equation 
(13).  If  one  simply  uses  a  double  summation  over  all  of  the 
near-field  data  for  each  K-space  point  the  computer  time 
quickly  becomes  exceedingly  large.  Another  less  time 
consuming  method  is  to  create  an  intermediate  step  in  the 
calculation.  This  requires  two  loops,  the  first  sums  over 
one  of  the  near-field  coordinates  for  each  value  of  the 
remaining  near-field  coordinate,  and  for  each  value  of  one 
of  the  K-space  coordinates.  The  second  loop  sums  over  the 
remaining  near-field  coordinate  for  each  K-space  coordinate 
pair  [15;  4:381]. 


S(ym»‘‘x)  =  E  E<Xn‘ym>  + 
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This  method  of  numerical  integration  is  analogous  to 
taking  the  one  dimensional  Fourier  Transform  of  the  x 
variable  for  each  value  of  y,  and  then  taking  the  one 
dimensional  Fourier  Transform  of  the  y  variable  for  each 
value  of  kj^  [14:118;  1:75;  7].  The  last  function  of  TRPROG 
is  to  load  the  calculated  vector  amplitude  values  into  a 
data  file  for  retrieval  by  the  DRPROG. 

DRPBQG 

This  program  calculates  the  first  partial  of  the 
electric  field  with  respect  to  z  by  implementing  equation 
(17).  The  values  of  the  vector  amplitude  function 
calculated  by  TRPROG  are  read  into  a  two  dimensional  array. 
Since  the  coordinate  transformations  were  performed  in 
TRPROG  this  is  a  considerably  easier  task.  Next,  the  double 
integral  of  equation  (17)  must  be  computed.  A  method  very 
similar  to  that  used  in  the  previous  program  is  used. 

Again,  an  intermediate  step  is  introduced  into  the 
computation.  The  first  loop  sums  over  all  values  of  one  of 
the  K-space  coordinates  for  each  value  of  its  real  space 
counter-part,  and  for  each  value  of  the  unused  K-space 
coordinate.  The  second  loop  sums  over  all  values  of  the 
remaining  K-space  coordinate  for  each  value  of  its  real 
space  counter-part,  and  for  each  value  of  the  first  K-space 


coordinate  used  [15]. 

S(kym,x')  =  ^ A(kjin»kyni) 
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The  major  difference  between  CSIPROG  and  TRPROG  is  that 
DRPROG  will  only  sum  over  the  vector  amplitude  values  of 
propagating  modes,  because  the  evanescent  modes  of  the 
hypothetical  test  antenna  are  assumed  negligible  at 
z  =  Zq  >  a  in  front  of  the  aperture. 

OQIVIPRQG 

The  combination  program  evaluates  equation  (19).  It  is 
an  extremely  simple  routine  which  is  separated  only  to 
accommodate  the  hardware  limitations  discussed  previously. 
Also,  to  reduce  the  need  for  large  storage  blocks,  the  near¬ 
field  data  originally  read  into  TRPROG,  is  re-read  into 
OOMPROG. 

QUTPRQG 

The  output  program  re-calculates  equation  (13).  At 
first  glance  it  would  appear  that  needless  programming  went 
into  the  creation  of  OUTPROG,  however  there  is  one  major 
difference  between  TRPROG,  and  OUTPROG.  TRPROG,  by 
necessity,  calculated  the  vector  amplitude  function  for  a 
grid  of  points,  equidistant  from  each  other,  in  K-space. 

This  format,  however,  is  not  conducive  to  presentation  as  a 
graph.  The  vector  amplitude,  or  far-field,  pattern  is 
usually  presented  as  a  function  of  the  spherical,  angular 


coordinates  and  <t>  •  Therefore,  OUTPROG  calculates 
the  vector  amplitude  function  for  points  equidistant  in 
spherical  coordinates.  Also,  there  is  no  need  to  calculate 
a  complete  grid  of  data,  as  a  single  cut  in  the  B  direction 
for  a  fixed  value  of  will  yield  all  of  the  information 
required  for  comparison.  This  is  due  to  the  fact  that  the 
far-field  pattern  of  the  test  antenna  is  circularly 
s  ymme  trie. 

OUTPROG  calculates  vector  amplitude  values  for 
4>  -  7r/4  and  theta  ranging  from  >30°  to  -1-30°  in  increments 
of  1/3°  ,  equation  (14)  gives  the  functional  relationship 
between  the  K  vector  and  0  and  4>  •  The  computation 
portion  of  OUTPROG  is  identical  to  TRPROG. 


To  simulate  errors  in  the  positioning  of  a  probe,  a 
slight  modification  is  made  in  the  data  read  by  TRPROG. 

Two  sets  of  near-field  values  are  calculated  by  the  database 
program.  One  is  a  slight  distance  in  front  of  the 
theoretical  plane  of  measurement  and  the  other  is  a  slight 
distance  behind  the  theoretical  plane  of  measurement.  As 
long  as  the  x  coordinate  of  a  storage  location  is  negative 
TRPROG  will  read  the  data  from  the  further  plane.  When  the  x 
coordinate  is  positive  TRPROG  will  read  the  data  from  the 
closer  plane,  (see  figure  4). 


s 
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This  is  a  convenient  model  for  positioning  errors  due 
to  its  ease  of  implementation  and  due  to  the  predictability 
of  the  effects  these  errors  will  have  on  the  far-field 
patterns . 

For  small  displacements  at  the  y-axis,  the  magnitude 
difference  between  the  two  halves  comprising  the  near-field 
data  will  be  negligible.  The  phase  difference  created  by 
propagation  will  not,  however,  be  negligible  (17:24).  This 
system  can  be  approximated  by  a  two  element  linear  array, 
and  the  theory  of  linear  arrays  leads  one  to  expect  a  shift 
in  the  position  of  the  main  beam.  Indeed,  this  shift  is 
observed  in  the  computation  of  the  far-field  patterns. 


V.  THE  DATA 

This  chapter  will  present  the  data  produced  by  the 
correction  routine.  Included  is  a  discussion  of  the 
validity  of  the  simulation  used  to  model  probe  positioning 
errors,  as  well  as  an  evaluation  of  the  performance  of  the 
transformational  programs,  TRPROG,  and  DRPROG. 

The  Model 

As  outlined  in  chapter  IV,  the  model  used  to  simulate 
erroneous  probe  positioning  consists  of  a  discrete  jump, 
equal  to  2Az,  in  the  measurement  plane  at  the  y-axis  (see 
figure  4).  It  can  be  argued  that  this  is  not  a  realistic 
simulation  of  the  positioning  errors  found  in  an  actual 
measurement  facility.  While  this  is  true,  the  errors 
represented  by  this  model  are,  most  pirobably,  more  severe 
than  those  of  a  real  near-field  system.  Thus,  if  these 
exaggerated  errors  do  not  present  undo  difficulties  to  the 
correction  procedure,  the  program  can  reasonably  be  expected 
to  perform  at  least  as  well  when  confronted  with  the, 
relatively,  benign  inaccuracies  of  a  high  precision  near¬ 
field  measurement  system. 

The  only  limitation  presented  by  this  particular 
correction  routine  is  that  the  mean  z  position  of  the  probe 
be  equal  to  the  theoretical  measurement  plane  z=Zq.  This 
ensures  a  proper  phase  reference  for  the  far-field  pattern. 


The  general  procedure  for  the  computation  of  a 
corrected  far-field  pattern  is  outlined  in  chapter  111, 
figure  2.  The  presentation  of  the  data  for  comparison, 
however,  requires  that  the  far-field  pattern  be  computed 
neglecting  any  positioning  errors.  This  is  accomplished  by 
supplying  the  program  OUTPROG  with  the  raw  data  before  it  is 
processed  by  the  correction  routine.  The  generation  of  a 
complete  set  of  data  that  facilitates  the  comparison  of  the 
corrected,  and  uncorrected  patterns  is  outlined  in  figure  5. 


Figure  5.  Computation  Procedure  for  Data  Presentation. 


The  computer  simulation  was  performed  for  an  antenna 
whose  diameter  is  two  meters  and  is  operating  at  2.4  GHz. 
Data  sets  were  generated  for  two  values  of  the  parameter 
z.  The  first  set,  shown  in  figures  6  and  7,  is  for  a 
positioning  error  of  /25  (S  mm).  Figure  6,  shown 
below,  is  a  comparison  plot  of  the  far-field  pattern 
computed  from  a  perfect  scan  (i.e.  no  positioning  errors), 
versus  the  plot  of  the  far-field  pattern  computed  from  a 
scanning  .plane  distorted  as  shown  in  figure  4.  The  plot  is 
of  a  cut  made  along  the  =  ■'/4  contour.  The  main  lobe 

has  been  shifted  by  approximately  1/3  degree. 


Figure  6.  Perfect  Scan  Pattern  vs.  Uncorrected  Pattern 


The  shift  in  the  main  lobe  closely  matches  the  shift 
predicted  for  a  two  element,  linear  array.  The  theoretical 
shift,  along  the  <t>  =  tt /4  cut  line  is  given  by: 
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Figure  7  is  a  plot,  made  over  the  same  contour,  of  the  far- 
field  pattern  of  a  perfect  scan  versus  the  far-field  pattern 
of  the  corrected  scan.  It  is  evident  that  the  main  lobe  has 
been  moved  back  towards  the  undistorted  position. 


CORRECTED 

PERFECT 


The  predicted  shift  in  the  main  lobe  is  given  by: 

^  ^  -i/47r/lo\  /n  \ 
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still  very  close  to  the  observed  shift. 

Figure  9  plots  the  perfect  pattern  versus  the  corrected 
pattern . 
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To  fully  understand  the  correction  procedure,  one  must 
observe  what  happens  to  the  near-field  magnitude,  and  phase 
patterns  of  the  corrected,  and  uncorrected  scans. 

Figures  10  and  11  are  plots  of  the  corrected  and 
uncorrected  phase  paterns  of  the  near-field  scans.  Figure 
10  is  a  cut  of  the  near-field  phase  distribution  as  the 
probe  moves  in  the  x  direction  on  the  distorted  scan  surface 
vs.  the  same  cut  made  on  the  perfect  scan.  Note  that  the 
jump  in  phase  across  the  origin  is  due  to  the  distorted 
scanning  plane. 
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Figure  10.  Near-Field  Phase,  Uncorrected  Scan. 


Figure  11.  Near-Field  Phase,  Corrected  Scan 
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In  figure  13  the  corrected  magnitude  pattern  is  shown, 
again  as  the  probe  moves  in  the  x  direction  along  the 
distorted  scan  surface  vs.  the  same  contour  on  the  perfect 
scan.  It  can  be  seen  that  the  pattern  differs  from  the 
uncorrected  scan  in  that  the  curve  appears  to  follow  the 
perfect  scan  slightly  more  closely  with  the  exception  of  a 
constant  shift  for  all  values  over  the  main  beam  of  the 
near-field  pattern.  Since  the  far-field  patterns  calculated 
from  this  data  were  normalized  with  respect  to  the  main 
beam,  this  constant  shift  does  not  seriously  affect  the  far- 
field  patterns. 


Figure  13.  Near-Field  Magnitude,  Corrected  Pattern. 
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There  are  two  main  components  of  the  transformation, 
and  correction  process,  the  program  TRPROG,  and  the  program 
ORPROG.  Both  routines  share  a  common  disadvantage  in  that  a 
complete  grid  of  data  must  be  computed  in  order  to  implement 
the  correction  routine.  This  is  due  to  the  fact  that,  for 
this  model,  positioning  error  corrections  must  be  performed 
for  each  point.  The  consequence  of  this  is  that  the 
computation  time  is  greatly  extended.  This  could  be  avoided 
by  using  a  more  efficient  transformation  algorithm,  such  as 
a  Fast  Fourier  Transform,  and  setting  some  tolerance 
threshold,  that  is  only  correct  those  measurements  that  are 
determined  to  be  excessively  inaccurate.  This  can  be  done 
using  these  techniques  since  the  correction  procedure  is 
point  by  point  in  nature.  In  this  way,  the  derivative  of  the 
near-field  values  would  only  have  to  be  computed  at  the 
points  where  it  is  needed  most. 

Another  point  to  consider,  is  the  accuracy  of  the 
presented  data.  The  far-field  patterns  were  computed  at 
increments  ot  6  =  1/3  degree,  then  normalized  to  the  largest 
value  computed.  In  this  way  some  error  is  introduced  into 
the  output.  If  the  actual  maximum  of  the  far-field  pattern 
occured  between  two  data  points,  the  program  would  tend  to 
mask  the  fact.  This  data  spacing  was  chosen  to  expedite  the 
execution  time  while  still  providing  reliable  accuracy.  If 
for  a  given  utilization,  the  main  lobe  was  all  that  needed 
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The  main  goals  of  this  investigation  were  accomplished 
in  that  the  effects  of  probe  positioning  errors  were 
investigated  and  found  to  be  significant,  even  for  errors 
much  less  than  a  wavelength.  The  viability  of  the 
correction  routine  was  investigated  and  determined  to 
adequately  correct  the  errors  introduced  by  inaccurate  probe 
positioning  given  that  those  errors  are  small  in  comparison 
to  a  wavelength. 

At  first  glance,  the  cases  studied  in  this  report  may 
seem  to  represent  an  extreme  case  of  inaccurate  probe 
positioning.  This  is  not  the  case  since  the  errors  were 
represented  with  respect  to  wavelength.  If  the  test  antenna 
had  been  operating  at  30  GHz  instead  of  2.4  GHz,  a 
positioning  error  of  A /2S  would  represent  a  distance  of 
0.4  mm,  a  very  fine  tolerance  when  it  has  to  be  maintained 
over  a  space  of  two  meters! 

The  use  of  this  correction  routine  could  reduce  the 
need  for  extremely  accurate,  and  therefore  expensive, 
positioners  without  sacrificing  reliable  far-field  patterns. 
Several  questions  remain  to  be  answered,  however,  before 
this  routine  could  be  implemented.  What  is  the  relationship 
between  the  correction  routine,  and  errors  in  the  magnitude 
distribution  of  the  corrected  near-field  scan?  Is  the 
correction  routine  responsible  for  the  magnitude  errors,  or 
are  they  due  to  the  crude  integration  algorithms  used? 


Can  the  routines  be  modified  so  that  they  can  be  used  in 
conjunction  with  scanning  geometries  other  than  planar? 

Along  with  these  technical  questionsi  there  are  several 
problems  involved  in  the  application,  and  optimization  of 
this  procedure.  The  programs  need  to  be  reworked  using  a 
more  efficient  Fourier  transformation  algorithm.  Probe 
correction  needs  to  be  taken  into  consideration.  The 
problem  of  the  presence  of  noise  in  a  real  data  set  needs 
to  be  investigated.  These  questions  represent  a  valuable 
source  of  research  opportunities  which  require  detailed 
investigations  and  could  lead  to  follow  on  theses  for  future 


AFIT  students 
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APPENDIX  A: 


The  following  is  a  listing  of  the  programs  used  for 
this  thesis.  They  are  presented  in  the  order  of  operation 
as  outlined  in  the  text. 


THE  DATABASE 

PROGRAM  DTPROG 

COMPUTES  THE  NEAR-FIELD  VALUES  FROM  THE  FAR-FIELD 
PATTERN  OF  A  PISTON  RADIATOR  USING  SPHERICAL 
HAMONIC  EXPANSION  (FOR  Z=CX>NSTANT) 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENS  ION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENS ION 


PN(IOOO) 

PNLOG(IOOO) 

PKROKZ(IOOO) 

PHASE(IOOO) 

XJKA02(500) 

XJMKA2(500) 

XJKR(IOOO) 

XNKR(IOOO) 

PL(2S0) 

BL(250) 

F(500) 

FF(500) 


INITIALIZATION 

PI=4.*ATAN(1. ) 
XLIM=1.2 

XKWA=(8.000001)*2.*PI 

ZOA=1.1375 

ZK=XKWA*ZOA 

MKR=3.*XKWA 

RNGEPN=-20. 

LMAX2=XLIM*XKWA/2. 

LMAX=LMAX2«2+9 

LMAXP1=LMAX+1 


CALL  SPHEBE((LMAX-t^l)/2,XKWA/2.  ,RES,XJKA02) 
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CALL  SPKEBE(-(LMAX+l)/2,XKWA/2. ,RES,XJMKA2) 
CALL  FACT(LMAX,XF,F) 

CALL  FCFACT((2*LMAX-1) ,XFF,FF) 


THIS  ROUTINE  COMPUTES  THE  SUMMATION  CONSTANT,  BL 

IF(LMAX.LT.l)  00  TO  100 

DO  20  L1=1,LMAX,2 

L=L1-1 

L02=L/2 

SQP=(FF(L+1)/F(L02+1)*XJKA02(1))**2*(-1)**L02 

IF(L-1)40,40,SO 

CONTINUE 

IF(L02.LT.l)  00  TO  105 

DO  30  IQP=1,L02 

F1=FF(L+1-2*1QP)/F(L02-IQP+1) 

F2=FF(L+1+2*1QP)/F(L02+IQP+1) 

ISUB=IQP+1 

F J=XJKA02 ( I  SUB ) *XJMKA2 ( I SUB ) 

QP=2.*Fl*F2*FJ*(-l)**(L02+IQP+2) 

SQP=SQP+QP 

CONTINUE 

CONTINUE 

CONTINUE 

BL(L1)=(2.*L+1)/2.**(L+1.)*SQP 

CONTINUE 

CONTINUE 


DO  200  IKRT=1,S01 
XKRT=MKR* (I KRT- 1 . ) / 5 0 0 
XKWR=SQRT ( XKRT* • 2  +ZK*  *  2 ) 

THETA=AOOS  (ZK/XKVm) 

CALL  SPHEBE(LMAXP1,XKWR,RES,XJKR) 
CALL  SPHERE (>LMAXP1,XKWR, RES, XNKR) 
IF(LMAXPl.LT.l)  00  TO  115 

DO  10  Isl,LMAXPl 

XNKRd  )=XNKR(I  +  1  )•(-!) ♦•  I 

CONTINUE 

CONTINUE 

CALL  LECEND(LMAXP1, THETA, PLL, PL) 
SUM1=0. 

SUM2=0. 

IF(LMAX.LT.1)00  TO  120 


c 


I 


DO  90  L1>1,LMAX,2 
L=L1-1 

TL11=L1*(ZK/XKWR*PL(L1)-PL(L1+1))*XJKR(L1) 
TL12=(L*XJKR(L1)-XKWR*XJKR(L1+1))*PL(L1) 
CSTHTA=ZK/XKWR 
TL12=TL12*CSTHTA 
SUM1=(TL11+TL12)*BL(L1)+SUM1 
TL21=L1*(ZK/XKWR*PL(L1)-PL(L1+1))*XNKR(L1) 
TL 2 2 = ( L • XNKR ( L 1 ) - XKWR • XNKR ( L 1 + 1 ) ) • P L ( L 1 ) 
TL22=TL22*CSTHTA 
SUM2=(TL21+TL22)*BL(L1)+SUM2 
90  CONTINUE 

C 

120  CONTINUE 

CX)EF= .  5  *XKWA*  •  2 /XKWR 
P=CI)EF  •  SQRT  ( SUMl  ♦  ♦  2  +SUM2  •  •  2 ) 

PN(IKRT)=P 

PNLOG ( I KRT ) = 2 0 . • ALOGl 0 ( PN ( I KRT ) ) 

PKROKZ( IKRT)=P*XKWR/ZK 
PHASE ( I KRT ) = ATAN2 ( SUM2 , SUMl ) 

200  CONTINUE 

C 

OPEN(UNIT=10,FILE=‘DAT2' ) 

C 

DO  18  1=1,501 
XI=(3./500,)*(I-1.) 

WRITE(10,17)  PHASE(I),PN(I) 

18  CONTINUE 

C 

17  FORMAT(E20.7,10X,E20.7) 

ENDFILE  10 
CLOSE(UNIT=10) 

STOP 

END 

SUBROUTINE  FACrr(N,F ,FMTRX) 

C  THIS  ROUTINE  COMPUTES  0  TO  N  FACTORIAL,  AND  STORES 
C  IN  THE  MATRIX,  FMTRX 

C 

DIMENSION  FMTRX(500) 

FMTRX(1)=1.0 

IF(N)11,11,20 

20  IF  (N.LT.l)  CO  TO  100 

C 

DO  10  1  =  1, N 

FMTRX( I  + 1 ) =FMTRX( I ) *  I 

10  CONTINUE 
C 

100  CONTINUE 

11  CONTINUE 
F=FMTRX(N+1) 

RETURN 


CONTINUE 

F=FMTRX(N+1) 

RETURN 

END 

SUBROUTINE  FCFACT(N,FF,FFMAT) 

COMPUTES  THE  FIRST  N  FAC-FACTORIAL  AND  STORES  IN  FFMAT 

DIMENSION  FFMAT(IOOO) 

FFMAT(1)=1.0 
IF(N)11,11,20 
IF(N.LT.1)00  TO  100 

DO  10  1=1, N, 2 

FFMAT( I +2 ) =FFMAT( I ) • I 

CONTINUE 

CONTINUE 

CONTINUE 

FF=FFMAT(N+2) 

RETURN 

END 

SUBROUTINE  SPHEBE(LEL,X, RESULT, T) 

COMPUTES  SPHERICAL  J  OF  ORDER  0  TO  LEL  AND  STORES  IN 
T 

DIMENSION  T(2500) 

IF(X)18,12,18 

IF(LEL)16,13,14 

RESULT=1.0 

T(1)=RESULT 

RETURN 

RESULT=0.0 

T(l)=1.0 

LELP1=LEL+1 

IF(LELP1.LT.2)  GO  TO  100 

DO  40  I=2,LELP1 

T(I)=0.0 

CONTINUE 

CONTINUE 

RETURN 

RESULT=1.0E+300*(-l)**(LEL+l) 

T{1)=1.0 

MI=-LEL+1 

IF(MI.LT.l)  GO  TO  105 

DO  30  1=2, MI 
T(I)=1.0E+300*(-l)**I 
CONTINUE 
CONTINUE 
RETURN 


1F(LEL)5S,51,19 

MO=LEL 

J0=2*IFIX(X) 

IF(MO-JO)2,21 ,21 
MO=JO 
M0=M0+11 
T(MO)=0. 

T(MO-l)=1.0E-285 

L0=M0-2 

F=2*(MD-1) 

M0=M0-3 

I2=M0 

F=F-2.0 

T(I2+1)=(F+1. )/X*T(I2+2)-T(I2+3) 

IF(I2)4,3,4 
12=12-1 
GO  TO  232 
F=SIN(X)/X/T(1) 

IF(M0.LT.l)OO  TO  110 

DO  5  J=l,MO 
T(J)=F*T(J) 

CONTINUE 

RESULT=T(LEL+1) 

RETURN 

RESULT=SIN(X)/X 

T(1)=RESULT 

RETURN 

LO=-LEL+l 

T(1)=SIN(X)/X 

T(2)=OOS(X)/X 

IF(LO.LT,3)GO  TO  115 

DO  6  J=3,LO 

T(J)=(-2.*FLOAT(J-2)+l)/X*T(J-l)-T(J-2) 

CONTINUE 

RESULT=T(LO) 

RETURN 

END 

SUBROUTINE  LEGEND (N, THETA, P.PMTRX) 

COMPUTES  LEGEND  POLY  OF  ORDER  0  TO  N  AND  STORES  IN  PMTRX 

DIMENSION  PMTRX(500) 

X=OOS (THETA) 

PMTRX(1)=1.0 

PMTRX(2)=X 

NP1=N+1 

IF(N-1)10,10,11 

CONTINUE 

IF(NP1.LT.3)00  TO  100 


20 

C 

100 

10 


DO  20  1=3, NPl 
XI  =  I 

PMTRXd )=(<2.*XI-3.)*X*PMTRX(I-1)-(XI>2. )• 
PMTRX(I-2))/(XI-l.) 

CONTINUE 


CONTINUE 


S'*. 


CONTINUE 

P=PMTRX(NP1) 

RETURN 

END 
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%  laERQG 


PROGRAM  TRPROG 

THIS  PROGRAM  IS  DESIGNED  TO  CALCULATE  THE  FAR-FIELDS 
ANTENNA  PATTERNS  OF  AN  ARBITRARY  TEST  ANTENNA,  FROM 
THE  NEAR-FIELD  VALUES  ON  A  SCANNING  PLANE  A 
SPECIFIED  DISTANCE  FROM  THE  ANTENNA 


REAL  PHSKSSO) 

REAL  PHS2(550) 

REAL  PNl(SSO) 

REAL  PN2(S50) 

REAL  KO 

COMPLEX  EX1(151,151) 
COMPLEX  Al(100,100) 
COMPLEX  SM(155,100) 
COMPLEX  XI ,XMI 
COMPLEX  CSUM 
COMPLEX  KT 
COMPLEX  KZ 
REAL  DELTAZ 
REAL  ZO 
REAL  LAMBDA 
REAL  KA 

LAMBDA=1. 0/8.0 
PI=4.*ATAN(1.) 

DELTAZ =.005 
ZO=1.125 
XI=(0. ,1.) 

XMI=(0. ,-l.) 

KO= 2 •PI /LAMBDA 

DK=(1.0/40.0)*KO 

DX=4/116.0 


THIS  LOOP  WILL  READ  IN  THE  NEAR-FIELD  DATA 


OPEN(UNIT=10,FILE='DAT1’ ) 
OPEN(UNIT=9,FILE=»DAT2' ) 
REWIND(UNIT=10) 

REWIND  (UNIT=9) 


DO  12  1=1,501 
READ(10,17)  PHSKD.PNKI) 

R£AD(9,17)  PHS2(1) ,PN2(1) 

CONTINUE 

FORMAT!  E20.7,10X,E20.7) 

CLOSE(UNIT=10) 

CLOSE ( UN I T= 9) 

DO  45  1=1,117 

DO  40  J=l,117 

JL=J-1 

IL=I-1 

R=SQRT((-2+JL*DX)**2.+(-2+IL*DX)**2.) 

IQ=INT(R*1000.)-10*INT(100.*R) 

TNR=500./3. 

IF(IQ.GE.5)IR=1NT(R*TNR+1) 

IF(IQ.LT.5)IR=INT(R*TNR) 

1R=IR+1 

IF(JL*DX-2.LE.O.O)GOTO  37 

ER=PN1(IR)*00S(PHS1(1R)) 

EI=PN1(IR)*S1N(PHS1(IR)) 

CO  TO  39 

ER=PN2(IR)*OOS(PHS2(IR)) 

EI=PN2( IR)*SIN(PHS2( IR) ) 
EX1(J,I)=CMPLX(ER,E1) 
OPEN(UNIT=10,FILE='NDAT’ ) 

WRITE! 10, 117)  REAL!EX1!J,I)),A1MAG!EX1!J,I)) 
FORMAT!  E20.7,1X,E20.7) 

CONTINUE 

CONTINUE 


THIS  LOOP  CALCULATES  THE  APPROXIMATE  VECTOR  AMPLITUDE 
USED  IN  THE  DDZEXl  ROUTINE 


DO  60  1=1,81 

DO  55  J=l,117 
CSUM= ! 0 . , 0 . ) 

DO  50  M=l,117 

KA=!-2+!M-l)*DX)*!-KO+!l-l)*DK) 
KT=CMPLX!0. ,KA) 

CSUM=CSUM+EX1 ! M , J ) ‘CEXP ! KT ) 

•DX**2. 


CONTINUE 


SM(J,I)=CSUM 

CONTINUE 

CONTINUE 

DO  75  1=1,81 

DO  70  J=l,81 
CSUM= ( 0 . , 0 . ) 

DO  65  M=l,117 

KA=(-KO+(J-l)*DK)**2.+(-KO+(I-l)*DK)**2. 
IF(KA.GT.KO**2)  GO  TO  65 
IF(KA.LE.KO**2)KZ=SQRT(KO**2-KA) 

KT=CMPLX( 0 . , ( - 2  +  (M- 1 ) •DX) • ( -KO+ ( I - 1 ) *DK ) ) 

KT=KT+KZ*ZO*XI 

CSUM=CSUM+ SM ( M , J ) * 

CEXP(KT)/(2.*PI ) 

CONTINUE 

A1(J,I)=CSUM 
OPEN(UNIT=7,FILE=»FRST* ) 

WRITE(7,150)  REAL(A1(J,I)),AIMAG(A1(J,I)) 
CONTINUE 

CONTINUE 


FORMAT(  E20.7,1X,E20.7) 


o  o  oonoo  oooon  oooo 
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PROGRAM  DRPROG 

THIS  PROGRAM  WILL  CALCULATE  THE  DIRIVATIVE  OF  THE 
TANGENTIAL  NEAR-FIELDS  WITH  RESPECT  TO  Z. 


COMPLEX  Al(100,100) 

COMPLEX  DDZEXl 

COMPLEX  SM(155,1S5) 

COMPLEX  CSUM 

COMPLEX  XMI 

COMPLEX  KX,KY,KZ,KT 

REAL  REA 

REAL  IMA 

REAL  ZO 

REAL  LAMBDA 

REAL  KA 

REAL  PI 

REAL  KO 

REAL  DK 

REAL  DX 


INITIALIZATION 


LAMBDAsl. 0/8.0 

PI=4.0*ATAN(1.0) 

ZO=1.125 

XMI=(0. ,-l .0) 

KO=2.0*P I /LAMBDA 

DK=(1.0/40.0)*KO 

DX=4/116.0 


THIS  LOOP  WILL  READ  IN  THE  APPROX.  VECTOR  AMP. 


OPEN(UNIT=7,FILE=»FRST' ) 
REWIND(UNIT=7) 

DO  12  1^1,81 

DO  11  J=1.81 
READ(7,17)  REA, IMA 
A1  ( J , I ) =CMPLX(REA, IMA) 
CONTINUE 


I 


CONTINUE 


FORMAT(  E20.7,1X,E20.7) 
CLOSE ( UN I T= 7) 


THIS  LOOP  WILL  CALCULATE  THE  ACTUAL  DERIVATIVE 


DO  30  1=1,117 

DO  25  J=l,81 
CSUM=(0. ,0.) 

DO  20  M=l,81 
ML=M- 1 
IL=I-1 
JL=J-1 

KA= (ML*DK-KO ) • • 2 . + ( JL*DK-KO ) • • 2 . 
IF(KA.GT.(.75)*KO**2.)  GOTO  20 
IF(KA.LE.( .75)*KO**2. )  KZ=SQRT(KO**2 .-KA) 
KX=CMPLX( 0 . , - ( I L*DX- 2 ) • (ML*DK-KO ) ) 
KT=KX+KZ*ZO*XMI 

CSUM=CSUM+A1 (M, J ) ‘XMI ‘KZ^CEXP (KT) / ( 2 . *P I ) 
CONTINUE 

SM(I,J)=CSUM 

CONTINUE 

CONTINUE 


DO  45  1=1,117 

DO  40  J=l,117 
CSUM= ( 0 . , 0 . ) 

DO  35  M=l,81 

ML=M-1 

ILsI-1 

KY=CMPLX( 0 . , - (ML*DK-KO) • ( I L*DX- 2 ) ) 
CSUM=CSUM+SM ( J , M ) *CEXP ( KY ) ’DK* • 2 . 
CONTINUE 

DDZEX1=CSUM 

OPEN(UNIT=9,FILE='THRD' ) 

WRITE(9,100)  REAL(DDZEXl) ,AIMAG(DDZEX1) 
CONTINUE 

CONTINUE 

FORMAT(  E20.7,1X,E20.7) 


PROGRAM  OOMBN 


THIS  PROGRAM  WILL  COMBINE  THE  ORIGINAL  NEAR-FIELD  DATA 
WITH  THE  CALCULATED  DERIVATIVE  OF  THE  NEAR-FIELD  DATA, 
AND  OBTAIN  AN  APPROXIMATION  OF  THE  NEAR-FIELDS  ON  THE 
THEORETICAL  PLANE  OF  MEASUREMENT. 


COMPLEX  EX1(1S5,155) 
COMPLEX  DDZ(155,15S) 
COMPLEX  EXO 
COMPLEX  SUM 
REAL  DER,DEI 
REAL  DELTAZ,DX 
REAL  PHSKSSO) 

REAL  PHS2(5S0) 

REAL  PNKSSO) 

REAL  PN2(550) 


DELTAZ=(1.0/8.0)*.l 

0X^4/116.0 


THIS  LOOP  WILL  READ  IN  THE  NEAR-FIELD  DATA 


OPEN(UNIT=10,FILE='DAT1' ) 
OPEN(UNIT=9,FILE=»DAT2’ ) 

R£WIND(UNIT=10) 

REWIND  (UNIT=9) 

DO  12  1=1,501 

READ(10,17)  PHS1(I),PN1(I) 

READ(9,17)  PHS2(1),PN2(1) 

CONTINUE 

FORMAT(  E20.7,10X,E20.7) 

CLOSE(UNIT=10) 

CLOSE (UNIT=9) 

DO  45  1=1,117 

DO  40  J=l,117 

JL=J-1 

IL=I-1 

R=SQRT((-2+JL*DX)**2.+(-2+lL*DX)**2.) 


IQ=INT(R*1000.)-10*INT(100.*R) 

TNR=500./3. 

IF(IQ.GE.5)IR=INT(R*TNR+1) 

IF( IQ.LT.5 ) IR=INT(R*TNR) 

IR=1R+1 

IF(JL*DX-2.LE.O.O)GOTO  37 

ER=PN1(IR)*00S(PHS1(IR)) 

E1=PN1(1R)*SIN(PHS1(IR)) 

CO  TO  39 

ER=PN2(IR)*CX)S(PHS2(IR)) 
EI=PN2(IR)*SIN(PHS2(IR)) 
EX1(J,I)=CMPLX(ER,EI) 
0PEN(UNIT=8,FILE=*NDAT’ ) 

WRITE(8,3S)  REAL(EX1(J,I)),AIMAG(EX1(J,I) ) 
CLOSE ( UN I T= 8) 

CONTINUE 

CONTINUE 


THIS  LOOP  WILL  READ  IN  THE  FIRST  DERIVATIVE, 


OPEN (UN I T=9 , F I LE= » THRD ’ ) 
REWIND(UNIT=9) 

DO  15  J=l,117 

DO  10  1=1,117 
READ(9,21)  DER,DEI 
DDZ ( I , J ) =CMPLX(DER , DE I ) 
CONTINUE 

CONTINUE 

CLOSE (UNIT=9) 

FORMAT(  E20.7,1X,E20.7) 


THIS  LOOP  WILL  COMBINE  EX1(X,Y)  AND  DDZ(X,Y) 


DO  25  J=l,117 

DO  20  1=1,117 
IL=I-1 

IF((IL*DX-2).LE.0.0)  GO  TO  18 
SUM=EX1 ( I , J ) -DELTAZ’DDZ (I , J ) 
GO  TO  19 

SUM=EX1  ( I ,  J )  -^DELTAZ^DDZ  ( 1 ,  J ) 
EXO=SUM 

OPEN(UNIT=7,FILE=»EXOUT» ) 


WR1TE(7,35)  R£AL(EXO) ,AIMAG(EXO) 
CONTINUE 

CONTINUE 

FORMAT(  E20.7,1X,E20.7) 

STOP 

END 


ouTPRnr. 


PROGRAM  OTPROG 


THIS  PROGRAM  WILL  CALCULATE  THE  PAR-FIELD  PATTERN 
OF  THE  TEST  ANTENNA  AND  PRESENT  THE  MATERIAL  IN  A 
FORM  SUITABLE  FOR  VIEWING. 


COMPLEX  EX0(1S5»15S) 

COMPLEX  AO(37S) 

COMPLEX  KA 
COMPLEX  CSUM 
COMPLEX  XI 
COMPLEX  K1,K2 
REAL  EFF(SOO) 

REAL  KX, KY,KZ,DX, PI, LAMBDA, KO, PHI ,DT,ZO 
REAL  ER(1S5, 155) ,£1(155,155) 

INITIALIZATION 


DX=4/116.0 

BIG^iO.O 

PI=4.0*ATAN(1.0) 

LAMBDA^l. 0/8.0 

KO=2.0*PI /LAMBDA 

ZO=1.125 

XI=(0.,1.0) 

PHI=PI/4.0 

DT=PI/540 


THIS  LOOP  WILL  READ  IN  THE  NEAR-FIELD  DATA 


OPEN(UNIT=7,FILE=»SCND») 

REWIND(UNIT=7) 

DO  12  1=1,117 

DO  11  J=l,117 

READ(7,17)  ER(J, I) ,EI(J,I) 

EXO(J,I)=CMPLX(ER(J,I),EI(J,I)) 

CONTINUE 

CONTINUE 

FORMAT(  E20.7,1X,E20.7) 

CLOSE (UNITS 7) 


•V 


CALCULATE  THE  OUTPUT  OF  THE  THESIS! ! ! ! I 


DO  30  M=l,181 
ML=M- 1 

CSUM=(0.,0.) 

DO  25  1=1,117 
IL=I-1 


DO  20  J=l,117 
JL=J-1 

KX=KO*OOS(PHI )*SIN(ML*DT- (PI /6 . ) ) 
KY=SIN(PHI)*(KX/OOS(PH1 )) 
KZ=K0*CX)S(ML*DT-(PI76.)  ) 
Kl=XI*KX*(JL*DX-2.) 

K2=XI*KY*(IL*DX-2.) 

KA=KZ*XI*ZO+Kl+K2 

CSUM=CSUM+EXD( J, I )*CEXP(KA)*DX**2 . /( 2 .‘PI ) 
CONTINUE 

CONTINUE 

AO(M)=CSUM 

£F1=CABS(A3(M)) 

EF 2 =CABS ( AO ( M ) ) *KX/KZ 

EFF(M)=SQRT(EF1**2+EF2**2) 

IF(B1G.LE.EFF(M))BIG=EFF(M) 

CONTINUE 


LETS  OUTPUT  SHALL  WE? 


DO  35  M=l,181 
ML=M- 1 

THETA=FLOAT(ML) • ( 1 . / 3 . ) - 30 . 
OPEN(UNIT=10,FILE=»LAST’) 
WRITE(10,100)  THETA, EFF(M) /BIG 
CONTINUE 

FORMAT(  E20.7,1X,E20.7) 
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In  recent  years,  near-fleld  antenna  measurement 
techniques  have  gained  a  good  deal  of  acceptance.  There 
are  several  errors  In  the  computed  faur-fleld  patterns  caused 
by  measiirement  Inaccuracies  In  the  near-fleld  data.  This 
paper  deals  with  deterministic  errors  Introduced  by  probe- 
posltlonlng  errors  In  planar,  near-fleld  scanning. 


By  utilizing  basic  near-fleld  theory,  as  well  as  a 
knowledge  of  the  positioning  errors.  It  Is  possible  to 
estimate  the  fields  at  the  correct  probe  position.  A 
computer  program  which  lessens  the  effect  of  probe 
positioning  errors  by  means  of  a  truncated  Taylor  series 
expansion  Is  used  to  demonstrate  this  improvement. 


To  simulate  inaccurate  probe  positioning,  a  database  of 
calculated  near-fleld  values  for  a  linearly  polarized, 
uniformly  Illuminated,  circular  aperture  was  ised.  Next, 
the  position  correction  program  was  applied  to  computer 
.generated  Inaccurate  data.  Finally,  far-fleld  patterns 
were  calculated,  and  compared  using  both  the  corrected  and 
uncorrected  data.  Results  of  the  comparisons  are  presented. 
Limitations,  and  areas  of  application  of  this  routine  are 
discussed. 
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